Chapter Two 


Transfer Function Approach 


In the previous chapter it has been indicated that modeling, analysis, and design 
of control systems can be performed in two domains, namely in the time and 
frequency domains. In this chapter we will consider the frequency (complex) 
domain technique, also known as the transfer function method. Our main goal is 
to present methods for finding the system transfer function. This is particularly 
important for systems composed of many blocks, where each block represents an 
internal transfer function. In Chapter 9, the frequency domain approach will be 
used to design controllers for linear time invariant systems. 

Modern control theory has its foundation in the state space approach; classical 
control theory is based on the transfer function approach. The state space 
method is widely used in modern control theory and practice due to the extensive 
support from modern packages for computer-aided control system analysis and 
design. The state space method will be considered in detail in Chapter 3. The 
transfer function approach is based on the Laplace and Z -transforms and their 
time derivative properties, which convert differential/difference equations into 
algebraic equations with complex coefficients. The algebraic equations obtained 
are frequency domain representations of the considered dynamic systems. The 
basics of the Laplace and Z -transforms are reviewed in Appendices A and B. 

In classical control theory, it is desirable to have a tool that permits analysis 
and design of control systems, especially when instead of knowing the internal 
state of the system, we just need to know the relationship between the system 
inputs and outputs. This can be facilitated if the system model is transformed 
from the time domain into the frequency domain. The transfer function—the main 
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concept of the frequency domain technique—is considered for both continuous- 
and discrete-time systems in Section 2.1. This book emphasizes continuous-time 
systems since many discrete-time results are easily derived by analogy from the 
corresponding continuous-time results. 

A conventional way of representing linear time invariant systems is via block 
diagrams. This provides a pictorial view of a control system. Block diagrams are 
considered in Section 2.2. Block diagram algebra is introduced in Section 2.3 as 
a suitable tool for obtaining transfer functions of systems whose block diagrams 
are known. The use of block diagram algebra to find the system transfer function 
is advisable for simple systems, but for complex systems it gets quite involved. 

The signal flow graph technique is employed in Section 2.4 as an alternative 
to the block diagram system representation. Mason’s gain rule, the main result of 
the signal flow graph technique, is an elegant way of finding transfer functions, 
especially for complex and high-dimensional systems. Several examples are 
given in order to demonstrate the power of Mason’s rule. 

In Section 2.5 we present specialized methods for finding transfer functions 
of sampled data control systems obtained by sampling continuous-time systems. 

At the end of this chapter, in Section 2.6, a laboratory MATLAB experiment 
on the system transfer function is designed. 

Chapter Study Guide and Objectives 

Students not completely familial - with the Laplace and Z -transforms should 
first read Appendices A and B. Instructors not interested in teaching transfer 
functions of sampled data control systems may skip Section 2.5 without loss of 
continuity. Sections 2.1-2.4 represent the core of the chapter. The main objective 
of this chapter is that students master a technique for finding transfer functions 
of any time invariant linear control system by using either the block diagram 
algebra or Mason’s rule. 


2.1 Frequency Domain Representations 

Real dynamic systems operate in real continuous time so that it is natural to 
describe and study their dynamical behavior and evolution in continuous time. 
This is done by using differential equations to model them. Some artificial 
dynamic systems operate in discrete time so that their models are represented 
by difference equations. In addition, discrete-time systems can be obtained by 
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discretizing continuous-time systems so that the obtained sampled data systems 
are also described by difference equations. The study of dynamic systems in both 
continuous and discrete time will be presented in detail in Chapter 3. 

Another way of studying continuous- and discrete-time systems is the fre¬ 
quency domain approach. This approach is performed in the space of complex 
numbers: by using the Laplace and 2-transforms, the differential/difference equa¬ 
tions are transformed into linear algebraic equations with complex coefficients. 
In general, it is easier to solve linear algebraic equations than linear differen¬ 
tial/difference equations, and hence the frequency domain approach seems to be 
very attractive. The frequency domain approach is often called the complex 
domain approach. Since all calculations have to be performed in the complex 
domain, and since the complex numbers s = a + ju and s = e~ sT (where 
T stands for a sampling period) are also known in engineering as complex fre¬ 
quencies, the common name for these methods is the frequency domain methods. 
The importance of such a representation of a system is especially emphasized in 
classical control system theory. 

We would like to point out that the frequency domain very often gives a better 
understanding of the actual control system phenomena than the time domain, but 
from the computational point of view the frequency domain is inferior to the time 
domain state space approach, especially for high-order dimensional systems. 


2.1.1 System Transfer Functions 

The system transfer function relates to the frequency domain system outputs and 
inputs. In other words, the system transfer function gives what is in between the 
system inputs and outputs, i.e. it indicates what kind of dynamic elements input 
signals have to face before they appeal - on the system outputs. This pictorial 
definition can be put in rigorous mathematical form by using the Laplace and 
Z -transforms. 

In the first part of this section we present the transfer functions for single¬ 
input single-output systems, and in the second part we study the general case of 
multi-input multi-output systems. 

Definition 2.1: The transfer function of a continuous-time single-input 
single-output system is defined as the ratio of the Laplace transform of the system 
output over the Laplace transform of the system input, when all initial conditions 


are zero. 
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Definition 2.2: The transfer function of a discrete-time single-input single¬ 
output system is defined as the ratio of the Z-transform of the system output over 
the Z-transform of the system input, when all initial conditions are zero. 

Note that there are several other ways to introduce the definition of the 
system transfer function (Franklin et al., 1990; Kuo, 1995), but all of them are 
basically the same. 

Some preliminary results on system transfer functions have been presented in 
Section 1.3. Since the presentation of the discrete-time transfer function parallels 
that for continuous time, we will mostly present the results for continuous¬ 
time transfer functions and give only the final results for discrete-time transfer 
functions. In Section 2.5 we will pay special attention to the transfer functions 
of sampled data systems. 

Consider a single-input single-output control system represented by an 
n -order differential equation, that is 


d n y(t) 
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( 2 . 1 ) 

with n > in. If the initial conditions are zero, its complex counterpart is obtained 
simply by a substitution of d l /df by s l and y(t) — Y(s), u(t ) — U(s) (see 
(a.4) in Appendix A), to give 


(s n -\- a n _i-s n 1 + • • • + (ii-s + (io)y (s) 

= {b-mS™ + b m -i-s m 1 + • • • + b\s + bo) U( s) 


( 2 . 2 ) 


Hence, the transfer function of this system is 

} ( s ) _ b m s m + b m -\s m 1 + • • • + b\s + bo _ N m {s) 

(S} ~U{T)~ s" + a n - lS "-i + ■ ■ ■ + a lS + ao ~ ~dJT) ( ’ 5) 

Similarly, for discrete-time systems the transfer function is obtained by applying 
the Z -transform to the difference equation describing system dynamics 

y(k + n) + (in-pyik + n - 1) +-f apyik + 1) + a 0 y(k) 

= b m u(k + m) + b m -\u(k + m — 1) + • • • + biu(k + 1) + bou(k) 


(2.4) 
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This yields the discrete-time transfer function of the form 

q, , _ 4 (t) _ b m z m + b m -\z m 1 + • • • + b\z + bo _ N m (z) 
U(z) z n + +-h a x z + (Iq D n (z) 


(2.5) 


Polynomials iV m and /!„ (with the s and s arguments dropped) have real 
coefficients, and for the so-called proper systems (real physical systems or causal 
systems) it must be satisfied that m < n. The meaning of a proper system is 
that the system cannot respond before an input to the system is applied (system 
causality, see for example, Kamen, 1990). 

The polynomial in the denominator of a single-input single-output system 
transfer function, I) „, is called the characteristic polynomial , and its roots are 
known as the system poles. At any of these n roots the denominator polynomial 
D n is zero, so that the overall transfer system function becomes infinite. The 
roots of the numerator polynomial N m are called the system zeros since at these 
m values both the numerator polynomial and the system transfer function are 
zero. If the system poles and zeros are known, the transfer function G(s) can 
be recorded in pole-zero form as 


G(s) = A 4s + Zl.)(.S+Z 2 )---(,S+Zm) 

{s + Pl)(s + p- 2 ) ■ ■ ■ (s + Pn) 

or in time constant form 

GW = t ,'Vtl)lvtl|'-(1.» + l) 

( T ai S + l)(r 02 S + 1) • ■■(T an S + 1) 


( 2 . 6 ) 


(2.7) 


From (2.3) and (2.6)-(2.7) we have 


K 


— bmi 


Kr 


-~1~2 

PlP'2 ■■■Pn 


T a , = —, i = 1,2,.^ n; r 6 _ = — , i = 1,2,...,???. 

Pt ~i 

Discrete-time transfer functions can be represented by the forms identical to (2.6) 
and (2.7) with the complex frequency -s replaced by the complex frequency ,r. 
Example 2.1: The discrete transfer function of the following system 


v(k + 4) + 3 y(k + 2) - y(k + 1) + 5 y(k) = u(k + 1) + 2 u(k) 



46 


TRANSFER FUNCTION APPROACH 


is obtained by assuming that all initial conditions are equal to zero and by applying 
the derivative (left shift in time) property of the Z -transform, which leads to 

(T 4 + 3 z 2 - z + 5)rU) = (z + 2)U(z) 


so that 


G(z) 


Y(:) 

U(z) 


z + 2 

z 4 + 3 z 2 — z + 5 


o 


Example 2.2: For the transfer function 

.s 3 + 0.4.s 2 - 0.95.S - 0.45 
■ ■ “ s 5 + 8.3s 4 + 23.1s 3 + 26.2s 2 + 10.4s 

we find the pole-zero form by using MATLAB function tf2zp. This produces 

(s — l)(.s + 0.5)(-s + 0.9) 

.s(.s + l)(.s + 1.3)(.s + 2)(s + 4) 


o 


Example 2.3: The transfer function for the linearized system from Example 
1.1 is easily obtained by setting all initial conditions to zero and taking the 
Laplace transform, which leads to 


G(s) 


0(fl 

U(s) 


-1 


.s 2 — 1 


o 


It is important to point out that the system transfer function carries informa¬ 
tion about the system impulse response. In general, we have 

Y(s) = G(s)U(s) 

Y(z) = G(z)U(z) ( ' j 


Using impulse delta functions as inputs (U(s) = l,U(z) = 1), we get 
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so that, in the time domain 

g(t) = = £ _ 1 {G'(.s)} 

g(k) = Z-\Y impulse (z)} = Z- 1 {G{z)} 


The impulse response is obtained simply by finding the inverse transformation 
(from the frequency domain to the time domain) of the corresponding system 
transfer function. 

For multi-input multi-output (multivariable) systems the definition of the 
transfer function is more general since vectors and matrices are involved in the 
system transfer function description. 

Definition 2.3: Transfer functions of multivariable systems relate the fre¬ 
quency representation of system vector inputs and system vector outputs assum¬ 
ing that all initial conditions are equal to zero, that is 


Y(.s) = G(.s)U(.s) 
Y (z) = G( -:)U( -:) 


( 2 . 11 ) 


Note that G(.s) and G(,r) are matrices of dimensions p X r since we have 
assumed that the number of inputs is r and the number of outputs is p, so 
that vectors U(.s),U(,r) are of dimensions r x 1 and vectors Y(.s),Y(,r) have 
dimensions p X 1. Due to this “matrix” nature, one has to be careful while 
relating inputs and outputs for multivariable systems. For single-input single¬ 
output systems, one can write 

Y(s) = G(s)U(s) = U(s)G(s) 

Y(z) = G(z)U(z) = U (z)G(z) 1 .' 


Flowever, for multivariable systems this commutativity does not hold since in 
that case we are dealing with vectors and matrices. 

We have seen in Section 1.3 that the transfer matrix of a system with r 
inputs and p outputs has the form 


G pXr (s) = 

~Gn{s) 

G- 2 i{s) 

Gu{s) . 
G 22(s) • 

G\ r (s)' 
G‘2t{ •- ) 


-Gpi(s) 

G p2 (s) • 

• df p 7 1 ( S ) _ 
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where the coefficients Gij denote the transfer functions between the j th input 
and the ith output, when all other inputs are zero, that is 

G ‘ jM = 1T77 u (2 ' 14) 

|all inputs except jth are set to zero 


The following example illustrates the procedure for finding the transfer function 
for an inverted pendulum, which can be viewed as a multivariable system with 
two outputs and one input. 

Example 2.4: Recall a linearized model of the inverted pendulum, given 
by (1.77), that is 


(mi + m 2 ) 


d 2 x(t) 
dt 2 


mil 


d 2 x(t) d 2 9(t) 


dt 2 


+ l- 


dt 2 


d 2 0(t) 
dt 2 

= gO{t) 


= m 


Let us first find the transfer function Q(s)/F(s). The above system of two 
equations can be simplified by eliminating d 2 x/dt 2 , which produces 


m 2 l 


d 2 0(t) 
dt 2 


(mi + m 2 )g0(t) 


— F(t) 


After taking the Laplace transform of this equation, the transfer function is 
obtained as 


Gn(s) 


0 (f) 

F(s) 


l 

m-A 



9_ 

l 


The second transfer function, X(s)/F(s), is obtained by taking the Laplace 
transform of the second equation, that is 


•s 2 A'(.s) + /.s 2 0(.s) = gO(s) 


so that 


X(,) = 


which implies 


Is 2 - g /-s 2 — , 

- 12 -“- 12 - 


r< , X _ ls2 ~ ( J r< , » 
V T 2ll- 5 ) —-o 11V 1 
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The same result could have been obtained by taking simultaneously the Laplace 
transform of both the equations 


!"(???! + m-2)s 2 mils 2 1 



-F(.s)- 

■ 

i 

CO 

CO 

__i 

_0(.s) _ 


0 


and then solving this system of algebraic equations with respect to A'(.s) and 
0(s), that is 




(mi + m-2 )s 2 

mils 2 

-l 

-F(.s)- 


Gu(s) 

tt(.s) _ 


s 2 

Is 2 - g _ 


0 


_G-2i(s) 


F(s) 
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Sometimes systems are so complex that playing with algebraic equations 
in the complex domain in order to obtain the system transfer function(s) is 
mathematically very involved. A graphical system representation in terms of 
either block diagrams or signal flow graphs will help us to develop systematic 
methods for finding the system transfer function(s). 


2.2 Block Diagrams 

A pictorial description is a very convenient way of representing dynamic systems. 
It gives a clear picture of all components of the control system and the flow of 
information (signals) in the system. Such a representation is called the system 
block diagram. In the following, we show how to use block diagrams in order to 
obtain information about input and output variables, the relationships between 
these variables, and how to get the transfer function(s). In some cases the 
block diagram is used just to represent the composition and interconnections 
of a system. 

The simplest possible block diagram of a single-input single-output system 
is represented in Figure 2.1, where U(s) and Y(s) are, respectively, the Laplace 
transforms of the input and output signals, and G(-s) is the block transfer function. 
This is consistent with the transfer function definition given in (2.3). The rule 
that connects (2.3) and the block diagram in Figure 2.1, and which is also valid 
for all block diagrams, can be formulated as follows. 

Main Block Diagram Rule: The output signal is a product of the transfer 
function of the given block and block’s input signal. 
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Figure 2.1: Block diagram of a general linear system defined in (2.3) 

The basic structure of a single-input single-output closed-loop system with a 
non-unit feedback is presented in Figure 2.2a, and with a unit feedback is given 
in Figure 2.2b. Applying the main block diagram rule and taking into account the 
directions of the flow of signals in the system as indicated by arrows, we have 

y(.s) = G(s)E(s) 

E(s) = U(s) - H(s)Y(s) 



Figure 2.2: Simple feedback structures 
Eliminating E(s) from the last equation, it follows 

y(.s) = G(s)U(s) - G(s)H(s)Y(s) 
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so that for a single-input single-output control system the closed-loop transfer 
function , denoted by M(s), is given by 


M(s) 


Y(s) 

U(s) 


G(s) 

1 + G(s)H(s) 


(2.15) 


Note that if the system loop is open, which is the case for H(s) = 0, we have 
the open-loop system transfer function 


5'W = OMEM = GMUis) => = GM ,2.16) 


A more general block diagram of a closed-loop control system is given in 
Figure 2.3. 



Figure 2.3: Basic structure of a control system 

The two principal components (blocks) in this diagram are the plant and the 
controller. The plant has its own dynamics represented by the transfer function 
G p (s). The puipose of the controller (often called the regulator) G r {s) is to 
“reshape” the dynamics of the plant, such that the overall system transfer function 
has the desired form. In Figure 2.3, the signal D(s) represents the potential 
system disturbance. The main roles of the controllers (regulators) are system 
stabilization, improvement of the system transient response, reduction of steady 
state errors, disturbance rejection, etc. These issues will be addressed in the 
subsequent chapters. Here, we study only the problem of obtaining the system 
transfer function from the block diagrams. 
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In this block diagram we have two inputs, U(s) and D(s), and one output, 
Y(s). Setting D(s) = 0, we get 

y(.s) = G p (s)Gr(s)E(s) 

E(s) = U(s) - Y(s) 

that is 


Y(.s) = G p (s)GAs)U(s) - Gp(s)Gr(s)Y(s) 

so that the system closed-loop transfer function is given by 

M(«) = 1 (/S) = G p {s)Gr{s) 

■ ■ U(S) 1 + Gp(s)G r (s) 


(2.17) 


It is interesting to find the transfer function from the system disturbance 
D(s) to the system output Y(s). By setting U(s) = 0, we get 


E(s) = —Y(s) 

Y(s) = G p (s)[D(s) + Gr(s)E(s)\ = G p (s)D(s) - G p {s)G r {s)Y{s) 

so that 


y(,s) = g p ( s ) 

D (.s) 1 -\- G p (-S ) G r (-S ) 


(2.18) 


Since it is not desirable for the disturbance to affect the system output, the 
magnitude of the corresponding transfer function 


G p (s) 

1 -|- G p (s)G r (s) 


(2.19) 


should be minimized as much as possible by a proper choice of the controller 

Gr(s). 

In the case of multivariable systems, the closed-loop transfer function can be 
found using the same technique, bearing in mind that the corresponding quantities 
are vectors and matrices. The closed-loop block diagram of a multivariable 
system is given in Figure 2.4. 
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Figure 2.4: Multivariable feedback system 

From this block diagram by using the main block diagram rule, to be applied 
strictly in the order: output = transfer function X input , it follows 

Y(,) = G(,)E(,) 

E(.s) = U(.s) - H(.s)Y(s) 

so that 

Y(.s) = G(.s)U(.s) - G(.s)H(.s)Y(s) 

which implies 

Y(.s) = [I + G(.s)H(.s)] _1 G(.s)U(.s) (2.20) 

The closed-loop multivariable control system transfer function relates the fre¬ 
quency representations of the system vector input and vector output. From 
(2.20) it is given by 

M(.s) = [I + Gi.siIIi.s:,] 'Gi.sj (2.21) 

For H(.s) = 0, we have the open-loop transfer function of a multivariable control 
system obtained from (2.20) as G(.s). 

In the next subsection we show how to perform modeling and construct a 
block diagram from mathematical equations describing system dynamics. We 
consider a model of a DC motor, which is frequently used in control systems. 
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2.2.1 Modeling and Block Diagrams of a DC Motor 

A DC motor is an electromechanical energy converter which converts electrical 
energy into mechanical energy. It is often used as an actuator in control systems. 
Figure 2.5 illustrates such a motor schematically. In this section we present the 
modeling of a DC motor and draw two block diagrams corresponding to two 
different working conditions. 



Figure 2.5: Schematic of a DC motor 

The basic equations of a DC motor (electric part) are obtained from 
Maxwell’s electromagnetic theory. The magnetic flux < f>(t ) is proportional to 
the field current ij{ I ) 

= hi flt) (2.22) 

The torque produced by a motor is proportional to the product of the magnetic 
flux <f>( /) and the armature current i a (/) 

T(t) = k- 2 <f>(t)i a (t) = /■• i/*_*C * I !/ 1 1 (2.23) 

The motor’s electromotive force (induced voltage), denoted by e m (t), is propor¬ 
tional to the product of the magnetic flux and the rotor shaft rotational speed u(t) 


Cm(f) — — kik^i flt)ij(t) 


(2.24) 
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The last three equations are valid if the values of <f>(t),T(t), and e m (t) are close 
to their nominal values. 

On the mechanical side, the torque T(t), developed by the motor, is balanced 
by the load and disturbance torques. T(t) is also related to the rotational speed 
jj(i) by the differential equation 

T(t) = Ti(t) + T d (t) = J^p- + F'oj(t) + T d (t) ( 2 . 25 ) 

where J is the combined load and armature mass moment, and F is the viscous 
friction coefficient. Ti(t) represents the load torque and T d (t) is a disturbance 
torque (T d is frequently negligible). 

Balancing the voltages in the field and armature windings, we obtain 

VfW = + RFfW ( 2 - 26 ) 

v a (t) = L a dla p + R a i a (t) + e m (t) 

1 (t) {121) 

= Lapp 1 + Raia(t) + k^ifdMt) 

The above set of equations is nonlinear due to the presence of the products 
if(t)i a (t) and if(t)u>(t). However, usually one of the currents is kept constant. 
For constant i /(/), we have the so-called armature-controlled DC motor; if i a (t) 
is constant, then the motor is said to be field-controlled. Mathematical models 
for these two regimes are different. They are presented below. 
Armature-Controlled DC Motor 

In this case if(t) = Ij 0 = const, so that v fit) = II I Balancing the 
voltages and torques, we obtain from ( 2 . 23 ), ( 2 . 25 ), and ( 2 . 27 ) 

v a(t) = L a — Pj— + L : i : ! / ! + k 4 u(t.) ( 2 . 28 ) 

i : 1 1 1 = •/ L- F'oj(t) + T d (t) ( 2 . 29 ) 

where 

— h\ha,I j 0 , — k d k'2l jo 
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Quantities r a = L a /R a and r m = .J/F are usually called the system time 
constants. The Laplace transform of the above system of equations produces 


Us)= (L aS + R,F- M kMS)] 

(2.30) 

n(s) = (,/. s + F) [ k ^ S ) T di S )] 

(2.31) 


The block diagram for this system is easily drawn by looking at equations (2.30) 
and (2.31) and using the main block diagram rule. This is shown in Figure 2.6. 



Figure 2.6: Block diagram for an armature-controlled DC motor 
Field-Controlled DC Motor 

In this case i a = I ao = const . the following differential equation is obtained 
from (2.23) and (2.25) 


k e if(t ) = J + F'uj(t) + T d (t), k 6 = hhFo (2.32) 

Taking the Laplace transforms of (2.26) and (2.32), we get 

Vf(s) = R f (T f s + 1 )If{s), r ^ (2 33) 

kclf{s ) = F(r m s + 1 )0( .s) + T d (s) 

where tj is a time constant. Rearranging these equations in the form 

1 


I f (s) = 


Rf(T f S + 1) 


Vf(s) 


n(s) = 


i 


F(r m s +1) 


[kelf(s) - T d (s)] 


(2.34) 
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and using the main block diagram rule, the block diagram is easily drawn and 
is represented in Figure 2.7. 


TJs) 


V f (s) ^ 

l/Rf 



T(s) T,(s) 

l/F 


XjS+ 1 


K 6 




Figure 2.7: Block diagram for a field-controlled DC motor 

Real control systems can have complex structures including several local 
feedback loops, and many inputs and outputs. Local feedback loops can arise 
for two reasons: they are either the result of the physical nature of the specific 
element in the system or the system as an entity, or they can be intentionally 
built in, with the aim of achieving a desirable performance for the system. No 
matter how complex the block diagram of a system is, it can be reduced to one of 
the basic structures given in Figures 2.1-2.3 by using the block diagram algebra 
rules. In the following we present the main results of the block diagram algebra. 


2.3 Block Diagram Algebra 

Block diagram algebra is a set of rules that facilitates modification and simplifi¬ 
cation of block diagrams. The rules of block diagram algebra for continuous-time 
systems are quite simple. They are based on simple principles of algebra that 
are used for writing input-output relations for the specific blocks in the block 
diagram. 

Cascade (serial) connection: The transfer function equivalent to a serial 
connection of n blocks with transfer functions G'i(-s), G^G), ...>G' n (.s), repre¬ 
sented in Figure 2.8, is given by 


G(s) = G'i(.s)G' 2 (.s) • --Gnis) = 1J Gi(s) 

8 = 1 


(2.35) 
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This rule is obtained by generalizing the main block diagram rule. 



Figure 2.8 A serial connection of n blocks 


Parallel (tandem) connection: The equivalent transfer function for such a 
connection representing a summation of signals, given in Figure 2.9, is obtained 
as 

n 

G(s ) = G'i( s) + G- 2 {s) + • • • + G n ( s ) = 'y " Gi(s) (2.36) 

8 = 1 



Figure 2.9: Parallel connection of n blocks 

Feedback connection: The simplest form of a feedback control system is 
given in Figure 2.2. For such a system connection the transfer function is given 
by (2.15). 
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Example 2.5: In this example we demonstrate the procedure for obtaining 
the transfer functions for cascade, parallel, and feedback connections by using 
MATLAB. Consider the transfer functions 


G\(s ) 


5 

s(s + 1 )(s + 2) 


G-2( s) 


s + 4 
s + 5 


We obtain the cascade (series) transfer function by the following sequence of 
MATLAB operators 


% define G1(s) 
zl= [inf;inf;inf]; 

% all three zeros of Gl(s) are at infinity 
P l=[0;—1;—2]; 
kl = 5; 

[nl,dl]=zp2tf(zl,pl,kl); 

% zp2tf maps zero-pole transfer function into 

numerator-denominator transfer function 
% define G2(s) 
n2= [ 1 4]; 
d2= [ 1 5]; 

% find the series connection 
[ns,ds]=series(nl,dl,n2,d2); 

% print 

printsys(ns,ds,'s') 

Execution of these operators produces the following result 

r i a \r< t a \ - __ 

2 " “ s 4 + 8.s 3 + 17s 2 + 10s 

The transfer function for parallel connection is obtained by using 


[np,dp]=parallel(nl,dl,n2,d2); 
printsys(np,dp,'s') 


This produces 


G'i(s) + G- 2 {s) 


s 4 + 7.s 3 + 14.s 2 + 13s + 25 
.s 4 + 8.s 3 + 17.s 2 + 10s 


The feedback connection is executed by 
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[nf,df]=feedback(nl,dl,n2,d2,-l); 
% -1 indicates negative feedback 
printsys(nf,df ,' s ') 

which leads to 


G'i(.s) _ 5s+ 25 

1 + G'i(s)G' 2 (s) ~~ .s 4 + 8.s 3 + 17s 2 + 15s + 20 


o 

In addition to algebraic formulas (2.35) and (2.36), block diagram algebra is 
complemented by several “geometric” rules. Two of them are given below. 

Moving Pick-Off Point: In some cases it is desirable to move a pick-off point 
in front or behind a block in a block diagram, such that the terminal signals do 
not change their values. Figure 2.10 shows the equivalent block diagrams for the 
cases before and after replacement of the pick-off points. 

Moving Summing Point: An adder or subtracter may be moved from one 
side of a block to another as Figure 2.11 illustrates. It is easy to show that 
the diagrams on the left-hand and right-hand sides are equivalent. We leave the 
proofs to the reader. 


X 



Y 

► 

Y 

-► 


X 



Figure 2.10: Moving pick-off point transformation 
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Figure 2.11: Moving adder/subtracter transformation 

Using formulas (2.35) and (2.36) and moving pick-off and summing point 
rules, additional block diagram algebra rules can be derived. A summary of some 
additional useful block diagram algebra “geometric” rules is given in Table 2.1. 

Next we solve two simple examples in order to demonstrate the procedure 
for finding the system transfer function from block diagrams. 

Example 2.6: The original block diagram, presented in Figure 2.12a, is first 
simplified by moving the adder in front of the block G'i (Figure 2.12b), then by 
interchanging adder and adder/subtracter (Figure 2.12c), and finally by finding 
the corresponding closed-loop transfer functions (Figure 2.12d). 


o 

Example 2.7: The block diagram from Figure 2.13a is redrawn in Figure 
2.13b in order to explicitly indicate block connections and signal flows. In the 
next step, presented in Figure 2.13c, the closed-loop transfer function of blocks 
G'i and H 2 is found and the pick-off point is moved in front of block GG Finally, 
two closed-loop transfer functions are found (Figure 2.13d) and their cascaded 
connection is evaluated (Figure 2.13e). 


o 
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Figure 2.13: Simplification of the block diagram in Example 2.7 







64 


TRANSFER FUNCTION APPROACH 


Using the rules for simplifying the block diagrams presented in formulas 
(2.35) and (2.36), Figures 2.10 and 2.11, and Table 2.1, and finding the corre¬ 
sponding transfer functions is relatively easy for simple systems. However, for 
complex systems the procedure can be quite involved since it requires drawing 
of many intermediate block diagrams before the final (simple feedback form) 
is reached. Example 2.8 demonstrates the required procedure for a complex 
feedback system. 

Example 2.8: The reduction of a complex block diagram for a system shown 
in Figure 2.14a is illustrated in Figures 2.14b—i. The process of reduction is 
pretty much self-explanatory from the corresponding Figures 2.14b—i. The above 
simplification is primarily done by using the established rules, but in addition, 
one has to use common sense, as was done in going from Figure 2.14b to Figure 
2.14c and from Figure 2.14c to Figure 2.14d (see Problems 2.12 and 2.13). 

The final expression for the transfer function is given by 

U(-s) = _ G1G3G4G2 + G 5 ) _ 

U(s) 1 + G- 2 (HiG\ + H2G3) + 473&'l&'3&'4(&'2 + G'5) — H1H2G1G2G3G5 

It can be seen from this particular - example that for complex systems the 
block diagram algebra produces the required answer after many redrawings of 
the original block diagram. 



Figure 2.14a: Block diagram of the control system for Example 2.8 
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Figure 2.14f—i: Simplification of the block diagram from Figure 2.14a 


o 

In the next section we present another method, based on signal flow 
graph theory, known as Mason’s rule, which is particularly efficient for high¬ 
dimensional and complex systems. 
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2.4 Signal Flow Graphs and Mason’s Rule 

Another way of finding transfer functions of linear time invariant systems rep¬ 
resented by their block diagrams is the so-called Mason’s rule, which is based 
on signal flow graph theory (Mason, 1953, 1956). Signal flow graph techniques 
are also used in several other areas of engineering and sciences (Robichaud et 
al, 1962; Rao and Koshy, 1991). 

There is an analogy between block diagrams and signal flow graphs. The 
main elements of the signal flow graph technique are nodes and branches, with 
the nodes being connected by branches. A branch is equivalent to a block in 
the block diagram and represents the transfer function between the nodes. A 
branch consists of input node, output node, and an arrow showing the signal flow 
direction. A transfer function is associated with each branch. A branch in a signal 
flow graph and its transfer function counterpart are represented in Figure 2.15. 


U G(s) Y 

o->-o 



Figure 2.15: Equivalent elements in a block diagram and a signal flow graph 

A node represents a signal. The basic rule for nodes is that a signal at a 
node is equal to the sum of signals coming into the node from branches. Note 
that signals leaving a node do not count. It is only important to pay attention 
to the signals coming into a node. A signal entering a node from a branch is 
equal to the signal from the input node of that branch multiplied by the branch 
transfer function. 

In Figure 2.16 a simple feedback block diagram and its signal flow graph 
are given. It can be seen from this figure that the expressions for two signals 
at two nodes are given by 


E = 1 X U - HY 
Y = GE 


(2.37) 
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which produces our familial - closed-loop result 

Y = GU - GEY =*► Y = —U (2.38) 

1 + GH 


U ~ E 


1 7 


G 


n* U J E 9 Y J Y 




° 

1_ 

H 

< - 

-H 


(a) 

(b) 


Figure 2.16: Equivalence between closed-loop systems 

In order to be able to present the method for finding transfer functions by 
the signal flow graph technique, we need to introduce the following terminology. 

Source node — a node at which signals flow only away from the node. 
Input signals are represented by such nodes. 

Sink node — a node at which signals flow only towards the node. These 
nodes represent output signals. It is customary to extract the inputs and outputs 
out of a signal flow graph by using additional branches whose transfer functions 
are equal to 1 (see Figure 2.16b). 

Path — a succession of branches from a source node (input) to a sink node 
(output) with all arrowheads in the same direction which does not pass any node 
more than once. The path gain is the product of all transfer functions in the path. 

Loop — a closed path of branches with all arrowheads in the same direction 
in which no node is encountered more than once. A source node cannot be a 
part of a loop since each node in the loop must have at least one branch into the 
node and at least one branch out of it. The loop gain is the product of transfer 
functions of the branches comprising the loop. 

Nontouching loops — two loops are nontouching if they have no common 
node. 
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The above notions are demonstrated on a signal flow graph presented in 
Figure 2.17. 



Figure 2.17: Example of a signal flow graph 

In this figure the input node is U and the output node is V . There are 
two paths connecting the input U and the output V . One of them has a path 
gain equal to G 1 G 2 G 3 G 4 G 5 ., and the other has a path gain of GxG&G^G^. 
There are four loops in this signal flow graph with gains — G4H4, — G3G4 H>. 
— G4G2G3G4G5H3, and —G4GQG4G5H3. Two loops, — G\H\ and — G3G4H2, 
do not touch each other, i.e. they represent nontouching loops. 

Note that the signal flow graph contains fewer elements than the correspond¬ 
ing block diagram. Signal flow graphs can be simplified by employing similar 
rules to those we have been using in block diagram algebra in order to deter¬ 
mine system transfer functions. However, we have seen from the examples in 
the previous section that the simplification procedure based on block diagram 
algebra is quite lengthy. 

In Mason (1953, 1956) an elegant and powerful formula for finding the 
transfer function between input and output nodes was derived. That formula is 
known as Mason’s gain formula and is given by 

1 N 1 

G(s) = — Pk^-k = ^-(-PiAi + P 2 A 2 + • • • + PjvAjv) (2.39) 
k =1 

where /'/, is the path gain for the /.th path, N stands for the number of paths, A 
is the determinant of the signal flow graph, and A k is the cofactor of path k. 
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A and A*, are computed as follows: 

A = 1 - (sum of all loop gains) + (sum of products of the loop gains of all 
possible combinations of nontouching loops taken two at a time) — (sum of products 
of the loop gains of all possible combinations of nontouching loops taken three 
at a time) + ■ ■ ■ . 

A k = value of A. for the part of the flow graph not touching the kth forward 
path. 

The application of formula (2.39) is demonstrated on the signal flow graph 
presented in Figure 2.17. We have already found that there are two paths with 
the corresponding gains 


Pi = GiG 2 G 3 G 4 G 5 , P 2 = GiGiGsGe 

and four loops whose loop gains are 

L\ = —GiHi, L - 2 = —G 3 G±H - 2 
L 3 = —G\G 2 G 3 GiGr,H 3 , L4 = —GiG^GhGqHz 

There are also two nontouching loops with gains L\ and L 2 . Then 

A = 1 — (Zi + L 2 + L 3 + L4) + L\L- 2 

Apparently, if we eliminate the path G 1 (J 2 (J:-(J\(Jr>, the remaining signal flow 
graph will have no loops left, so that = 1. The same conclusion is obtained if 
we eliminate the path G 1 Gr.G iGr, leading to A 2 = 1. Thus, the transfer function 
of the considered signal flow graph, according to formula (2.39), is given by 

y(,s) _ C( . __ Pi x 1 + p 2 x 1 __ 

U{s) J 1 — (Zi + L 2 + L 3 + L 4 ) + L\L 2 


G'iG' 2 G' 3 G' 4 G'5 + G1G4G5G6 


1 + G\H\ T G 3 G4H- 2 T G\G 2 G 3 G4G 3 H 3 T G 1 G 4 G. 5 G 6 P 3 T G\G 3 G4H\H- 2 


In the above expression Gi,i = 1,...,6, and Hj, j = 1,2,3, are either constants 
or functions of the complex frequency -s. 

Note that all signals in a signal flow graph are mutually related by linear 
algebraic equations. For example, in the case of the signal flow graph given in 
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Figure 2.17 we have 

E 1 = U- H\E 2 - H 3 Y 
E 2 = G\E\ 

E 3 = G2E2 - H 2 E 5 

e 4 = g 6 e 2 + g 3 e 3 

e 5 = G4E4 
Y =G 5 E 5 

where Gi,i = and Hj,j = 1,2,3, are coefficients. By playing simple 

algebra with the above system of linear equations one is able to obtain the required 
relationship between Y and U, i.e. the required transfer function. However, that 
approach is not systematic. Mason’s formula (2.39) is derived by using Kramer’s 
determinant method for solving systems of linear algebraic equations, which are 
obtained by relating signals in a signal flow graph. The complete proof of formula 
(2.39) is beyond the scope of this textbook; it can be found in Mason (1956) and 
Younger (1963). 

Next, we give an example to find the system transfer function by using both 
block diagram algebra and Mason’s rule. 

Example 2.9: Consider the block diagram given in Figure 2.18. 



Figure 2.18: Block diagram of a feedback control system 


The block transfer functions are given by 

.si.sm,- g ' 3( - s)=2 
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iTi(.s) 


s 

5 + 4’ 


H 2 (s) 


5 s 

s + 2 


(a) Using block diagram algebra rules, this block diagram is simplified as 
shown in Figure 2.19. The required transfer function is given in Figure 2.19d. 



(a) 


(b) 


(c) 


(d) 


Figure 2.19: Block diagram simplification using block diagram algebra rules 


(b) The signal flow graph for this example is represented in Figure 2.20. It 
contains two paths and three loops. In addition the loops (J \ 11 \ and (J:- /!> are 
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nontouching ones. Expressions for path and loop gains are 
Pi = G\G-2f P2 = G3, L\ = — G'l-ffi, L2 = —G3H2, L3 = —G1G2H2 
The graph’s determinant and path’s cofactors are obtained as 

A = 1 — (Zi + L2 + X3) + L1L2, Ai I, A2 = 1 — Xi 


G 3 



Figure 2.20: Signal flow graph for the system given in Figure 2.18 

The required closed-loop transfer function, according to Mason’s formula (2.39), 
is given by 

p(-y) = Pi + P 2 (i-x 1 ) 

U{s) 1 — (L\ + L 2 + L 3 ) + L 1 L 2 

_ G1G2 + G3 + G1G3H1 

~ 1 + G\H\ + G3H2 + G'iG' 2 P 2 + G1G3P1P2 


After substitution of the given values for Gi(s),i = 1,2,3, and Hj(s),j = 1,2, 
the transfer function is obtained as a ratio of two polynomials with respect to 
the complex frequency -s. The final result for the transfer function is given in 
Part (c) of this example. 

(c) The same problem can be solved by MATLAB using its functions for 
feedback, parallel, and series connections. Note that the procedure given below 
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cannot be applied to any signal flow graph. It can be applied only to those 
with explicitly distinguished feedback loops, series, and parallel connections. 
However, using the SIMULINK package one is able to obtain the transfer function 
for any block diagram and the corresponding signal flow graph. The transfer 
function of the feedback control system given in Figure 2.18 is found by using 
the following sequence of MATLAB instructions 

% feedback configuration of G1 and HI 

[n,d]=feedback([0 0 5],[1 1 0],[1 0], [ 1 4],-1); 

% cascade connection to G2 
[n,d]=series(n,d, [0 2 ] , [ 1 0 ] ) ; 

% parallel connection to G3 
[n,d]=parallel(n,d,[2],[1]); 

% feedback connection with H2 
[n,d]=feedback(n,d,[5 0],[1 2],-1); 
printsys(n,d,'s') 

This MATLAB program produces the following result 

Vi.s! _ 2 .s' + I l.s' + 38. s 3 + Ki.s 2 + 60s + 80 
U(s) ~ 11s 5 + 57s 4 + 109s 3 + 68s 2 + 200s 

o 


Example 2.10: Consider the block diagram given in Figure 2.21. 



Figure 2.21: Block diagram of a feedback control system 
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The corresponding signal flow graph is presented in Figure 2.22. In this example 
we have two paths and five loops. There are no nontouching loops. 



Figure 2.22: Signal flow graph for the system given in Figure 2.21 

The corresponding path gains Pi,i = 1,2, loop gains Lj,j = 1,2,..,*5, signal 
flow graph determinant, and cofactors arc given by 

p _ r 1 r 1 r' r 1 p _ r 1 r 1 r 1 r 1 

ir\ — Cr i (jr 2 Or 3 (jr 4 ? ^2 — 

L\ = — G4G2H4, L'2 = — G2G3H2, Ls = G1G5G3H2G2H1 

L4 = — G4G5G3G4H3, X.5 = —G4G2G3G4B-3 

A 1 = l, A 2 = 1 

A = 1 + Gi G 2 -Hi + G2G3 H 2 + G'i G 2 G3G4H3 -\- Gi G.5 G4H3 — G\ G.5 G3H2G2H4 


so that the closed-loop transfer function for this system is obtained as 

YW = 

U(s) 


G1G2G3G4 + G1G3G4G5 


1 + GiG r >H\ + G2G3H2 + G1G2G3G4H3 + G1G3G4G5H3 


G1G2G3 G5H1H2 


o 
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Note that the same problem is studied in Example 2.8 using block diagram 
algebra. Comparing the required calculations done in Examples 2.8 and 2.10, it 
is obvious that for complex systems Mason’s rule is much more efficient than 
the block diagram algebra approach. 

It should be pointed out that Mason’s rule is also applicable to the signal flow 
graphs corresponding to multi-input multi-output linear systems (see Problem 
2.17d). Mason’s rule can be also used for finding transfer functions from any 
two nodes in the signal flow graph (see Problem 2.14). 

Having obtained the system transfer functions, we will be able to design 
controllers in the frequency domain such that feedback systems satisfy certain 
specifications, like desired transient and steady state responses. A frequency 
domain controller design technique based on Bode diagrams will be presented in 
Chapter 9. Note that Bode diagrams, in fact, represent the frequency plots, for 
•s = j'oj, of the magnitude and phase of the system transfer function. 

For discrete-time systems that are inherently discrete, duality can be em¬ 
ployed and the same rules for finding discrete transfer functions as for continuous¬ 
time systems are valid. However, there are some differences in the case of 
discrete-time systems obtained through sampling (sampled data systems). Trans¬ 
fer functions of sampled data systems are considered in the next section. 

2.5 Sampled Data Control Systems 1 

In determining discrete transfer functions of sampled data systems the procedure 
is a little bit more complex. In some cases, the corresponding transfer function 
even does not exist since it is impossible to find a linear relationship in the 
frequency domain between output and input signals. 

In Sections 2.5.1 and 2.5.2 we present procedures for finding the basic 
transfer functions of open-loop and closed-loop sampled data control systems. 
Section 2.5.3 studies the closed-loop transfer function for a special class of 
sampled data control systems known as digital computer controlled systems. 
Here, we present only the basics. For more information about the sampled data 
control systems, the reader is referred to specialized books (e.g. Astrom and 
Wittenmark, 1990; Ogata, 1987; Franklin et al, 1990; Kuo, 1992; Phillips and 
Nagle, 1995). 

1 This section may be skipped without loss of continuity. 
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2.5.1 Open-Loop Transfer Functions 

Depending on where the samplers are positioned, even for very simple block 
diagrams, several interesting cases may arise. These arc illustrated in Figure 
2.23, where T stands for the sampling period. 



Figure 2.23: Possible cascade connections in a sampled data system 


From Figure 2.23a, it follows that 

Y(s) = G(s)U*(s) 

Y*(s) = [G(s)F*(s)]* = G*(s)U*(s) 


(2.40) 


Formula (2.40) indicates one of the main properties of the starred Laplace 
transform defined in Appendix B in (b. 11), see also (b. 16). Since the starred 
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Laplace transform of a signal is equal to its 2-transform (see (b.12) in Appendix 
B), it follows from the above equation that 


y(~) = G(z)U(z) 


YU) 

U(z) 


G(z) = G*(s) ls=Ll ns 

I t 


(2.41) 


where T is the sampling period, and G( z ) the required transfer function. 

Serial connection of two cascaded blocks with samplers, presented in Figure 
2.23b, can be described by the following set of equations 


Ui(s) = Gi(s)U*(s) 
Y(s) = G 2 {s)Uf(s) 


(2.42) 


Using the same procedure as in the previous example, we have 

U*(s) = G${s)U*(s) 

Y*(s) = G* 2 (s)iq(s) 

and 

lh(z) = G 1 (z)U(z) 

Y(z) = G 2 {z)lh{z) 

which yields 

Y{z) = G 2 (z)Gi(z)U(z) 


so that the transfer function for this system is 

y = G 2 (z:)Gi(z) 


(2.43) 


(2.44) 


(2.45) 


(2.46) 


Note that in this case the transfer function is equal to the product of the transfer 
functions of each block, as in the case of two cascaded blocks representing 
continuous-time systems. 

A sampled data system with cascaded elements and no sampler in between 
is given in Figure 2.23c. Here, we have 


Y{s) - Gi{s)G 2 {s)U*{s) 


( 2 . 47 ) 
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The starred Laplace transform of (2.47) gives 


y*(s) = G\G* 2 (s)u* (s) 


(2.48) 


where G'lG'^-s) stands for [Gi(s)G 2 (s)]*. Then, it follows 


y(~) = GiG 2 (z)U(z) 
Y(z) 


U(z) 


= G'iG 2 (~) 


(2.49) 


It is important to note that the transfer functions in (2.46) and (2.49) are not the 
same, i.e. in general 


Gi(z)G- 2 (z) f G 1 G 2 G) (2.50) 

The last case, given in Figure 2.23d, is a serial connection of two elements 
with a sampler in between. For such a structure, we have 

U 1 (s) = G 1 (s)U(s) =* Ui(s) = G\U*(s) 

Vl.-O G 2 l.sj/ , l.s! 

Equation (2.51) gives, after Y(s) is staiTed 

y*(s) = G 2 (-s)G'iG*(.s) (2.52) 

so that 

y(s) = G2(z)Gilh(z) (2.53) 

It can be seen from (2.53) that, in this case, we are not able to identify the 
quantity that relates the system input and output in the frequency domain, in 
other words, for this particular open-loop sampled data structure the transfer 
function does not exist. 

The above discussion suggests that similar rules for continuous- and discrete¬ 
time block diagrams of the same structure are valid if there exists a sampler in 
front of each block of a sampled-data system (see Figure 2.23b). 
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2.5.2 Closed-Loop Transfer Functions 


Typical structures for closed-loop sampled data control systems are given in 
Figure 2.24. We will find that in two cases the closed-loop transfer function 



Figure 2.24: Four possible positions of a sampler in a closed-loop system 
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does not exist, and that in one case it has the same form as for the closed- 
loop continuous-time system given by (2.15). However, in this case, the transfer 
function depends on the complex frequency ' . 

According to Figure 2.24a, we have 



Y(.s) = G(s)E*(s) 

E(s) = U(s) - Yi(.s) = U(s) - G(s)H(s)E*(s) 

(2.54) 

By taking the starred Laplace transform, we obtain 



E*(s) = U*(s) - GH*(s)E*(s) 

(2.55) 

or 

1 + GH*(s) 

(2.56) 

Since 

Y*(s) = G*(s)E*(s) 

(2.57) 

it follows 

r 'W=l twt / 1 * 1 

(2.58) 


In terms of the Z -transform notation, the output is given by 


= TTmm u ^ (2 - 59) 


and the transfer function for the closed-loop system in Figure 2.24a is 


Y(G _ G(S) 

U(z) 1 + GH(z) 


(2.60) 


In a similar way, we obtain for the system in Figure 2.24b 


Y(z) 


GU(z) 

1 + GH(z) 


(2.61) 


which means that it is not possible to detennine the transfer function for this 
sampled data feedback system configuration. The same expression for the output 
as the one in (2.61) is obtained for the system given in Figure 2.24c. 
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Using the above procedure it is easy to find that the closed-loop system in 
Figure 2.24d has the transfer function given by 

Y{s) = G(z) 

U(z) 1 + G(z)H(z) K ’ 

which corresponds to (2.15), obtained for continuous-time closed-loop systems. 

The study of this section shows that the position of the sampler has a very 
important role, because it determines whether or not the input signal U(s) 
can be separated from the system dynamics. Therefore, we have seen from 
these examples that the transfer function exists if the sampler is just behind the 
comparator (subtracter), and conversely, does not exist if the sampler is in any 
other place in the system. 

It is important to note that in the case of sampled data systems, it is forbidden 
for a sampler and a continuous-time system element in cascade to mutually 
interchange their positions, because in that case the performance of the control 
system is changed. In other words, the commutative law is not, in general, 
applicable for such systems. 


2.5.3 Transfer Functions of Digital Control Systems 

A typical block diagram of a digital control system (digital computer controlled) 
is shown in Figure 2.25. One part of this configuration is a plant representing 
a continuous-time process that has to be controlled. Another part is a digital 



Figure 2.25: Block diagram of a digital control system 
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controller, usually a digital computer. The interfacing of these two parts is 
realized by a sampler (A/D converter), which converts analog signals into digital 
ones, and a zero-order hold (D/A converter), which converts digital signals into 
continuous-time signals. The continuous-time signals obtained are fed to the plant 
as control inputs. The block diagram showing transfer functions of all blocks in 
the above digital control system is given in Figure 2.26. 



Figure 2.26: Transfer functions in a digital control system 

In this system the error signal e( I ) = u( t ) — y( t ) is sampled so that the analog 
signal e(t) is converted into a digital signal by an A/D device. The digital form 
of the error signal e(k) is fed to the digital controller, whose transfer function 
is Gd.(z). After the controller has solved the difference equation described by 
this transfer function, the control signal v(k) is fed to a zero-order hold. The 
zero-order hold in time domain is represented by a unit pulse of duration T so 
that its the transfer function is 

G h (s) = £{h(t) - h(t - T)} = i - ~ e ~ Ts = - e~ sT ) (2.63) 

The plant transfer function is G p (s). From Figure 2.26 it follows 

Y{s) = G h (s )G P ( s )G* d ( s )E*(s) = G( s )G* d ( s )E*(s ), G( s) = G h (s )G p (s) 

(2.64) 

By the property of the starred Laplace transform given in (2.40), we have 

y*(.s) = G*(s)G* d (s)E*(s) (2.65) 
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which in the s-domain produces 


y ( ~) = G(z)G d (z)E(z) 


( 2 . 66 ) 


Since 

E(z) = U(z) - Y(z) (2.67) 

the last two equations give the required closed-loop transfer function of the system 
in Figure 2.26 as 


Y(~) _ G d (z)G(z) 
U(z) 1 + Gd(z)G(z) 


( 2 . 68 ) 


Note that in this introductory control theory course we will not pay attention 
to sampled data control systems. Material presented in this section is used to 
demonstrate the procedures for finding transfer functions of a class of linear 
discrete-time control systems known as sampled-data control systems and to 
indicate that special care has to be taken while finding the corresponding transfer 
functions. However, in this book we will study fundamental concepts and 
methods for discrete-time linear control systems mostly by using dualities with 
continuous-time linear control systems. 


2.6 Transfer Function MATLAB Laboratory Experiment 

Part 1. Consider the continuous-time system represented by its transfer 
function 

G(s)= , 3 + 1 

.s 2 + 5 s + 6 

Using MATLAB plot: 

(a) The impulse response of the system. Use impulse (num, den, t) with 
t = 0 : 0.1: 5. 

(b) The step response of the system. Use the function step (num, den, t ) 
with t = 0 : 0.1 : 5. 

(c) The system output response due to the input sin (2/). Use the function 
lsim (num, den, u, t ) with t = 0 : 0.2 : 2 0 and u=sin (2*t). 

(d) The system output response due to the input e~ t . Use the function 

lsim (num, den, u, t ) with t=0 : 0.1: 5 and u=exp (-t). 
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Part 2. A discrete-time system is described by the transfer function 


G(z) 


z — 2 

z 2 - 2.5s + 1 


Using MATLAB plot: 

(a) The impulse response of the system. Use dimpulse (num, den) and 
axis ( [0 10 0 1.5]). 

(b) The step response of the system. Use dstep (num, den) and axis ( [0 
10 0 3]). 

(c) The response due to the input si n (2/,- ). Use dlsim (num, den, u) with 
k=0:1:50 and u=sin(2 *k). 

(d) The steady state response due to the input 2~ k . Use dlsim (num, den, u) 
with u=2 . ~ (-k) and axis ([0 11 0 1.1]). Note that after 2 in 
MATLAB indicates a pointwise operation. 

Part 3. Consider a flexible beam system (Qiu and Davison, 1993) whose 
linearized model has the transfer function 


1.65s 4 - 0.331s 3 - 576s 2 + 90.6s + 19080 
■ ■ “ s 6 + 0.996s 5 + 463s 4 + 97.8s 3 + 12131s 2 + 8.11s 


Use MATLAB in order to find: 


(a) The system open-loop poles and zeros. Use the functions roots (num) 
and roots(den). 

(b) The system closed-loop transfer function assuming unity negative feedback. 
Use [numc,dene]=cloop(num,den,- 1 ). 

(c) The closed-loop poles and zeros and compare them to the open-loop poles 
and zeros found in (a). Use roots (numc) and roots (dene) . 

Part 4. The block diagram of a simple positioning control system using a 
field-controlled DC motor is shown in Figure 2.27. 

(a) Using the rules of block diagram algebra reduce this system to the basic 
feedback system shown in Figure 2.2a, and find the system transfer function 

Y(s)/U(s). 

(b) Find the transfer function of this system using Mason’s gain formula. 

(c) If K = 2, k\ = —0.05j k 2 = 0.16, and k 3 = 0.24, write the MATLAB 
script to find the transfer function of this system (see Example 2.5). 
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(d) Using MATLAB functions, find the poles of the closed-loop system. 

(e) For the given system poles find the partial fraction expansion using the 
MATLAB function residue. From the partial fraction expansion find 
analytically the system response to a unit step input. Check the obtained 
results using the MATLAB function step. 

(f) Plot the unit step response of the system for the time interval of 10 seconds. 



Figure 2.27: Position control system 
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2.8 Problems 

2.1 A control system is described by the following set of equations 

= <?(./), e ( f ) = u(t) - y(t) 

(a) Find the transfer function of this system and its impulse response. 

(b) Using the inverse Laplace transform, find the response of the system 
to a unit step input and zero initial conditions. 

(c) Using the final value theorem, calculate the steady-state error due to 
a unit step input, i.e. find e(t) for t — o o. 

(d) Plot the unit step response of this system by using the MATLAB 
function step (num, den, t) for the time interval of 10 seconds. 
Take t=0 : 0.1:10. 

2.2 Consider the electrical circuit given in Figure 2.28. 

(a) Find the voltage transfer function Vo(s)/Vi(s). 

(b) Suppose that an inductor Li is connected in parallel to resistor !( >■ 
Find the voltage transfer function for the modified circuit. 
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(c) If a constant input voltage r,( /) = 5 V is applied to the circuit, find 
the output voltages r,,( I j for the circuits given in (a) and (b) by using 
the inverse Laplace transform. 

(d) Using the final value theorem of the Laplace transform, find the steady 
state values of the output voltages in cases (a) and (b). 

(e) If R\ = l,C'i = l,i ?2 = 2,Zi = 1, plot the output v 0 (t) for the 
circuits given in (a) and (b) by using MATLAB. 


Ri 



Figure 2.28: An RC network 

2.3 Find the transfer function of the following continuous-time system 

y( 5 ) + 3 + 2 + y^ + 5y^ + y = 3 + u 

2.4 Consider the continuous-time linear system represented by the transfer 
function 

r = _ 2.s 5 + .s 3 - 3-s 2 + s + 4 _ 

■ 5.s 8 + 2s T - s 6 - 3.s 5 + 5.s 4 + 2.s 3 - 4.s 2 + 2.s - 1 

Use MATLAB to find: 

(a) The zeros and poles of the system. 

(b) The inverse Laplace transform of the transfer function. Use the 
residue function. 

(c) The system closed-loop transfer function assuming a negative unit 
feedback and find the corresponding closed-loop poles. 
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2.5 Consider the control system given in Figure 2.29. 



Figure 2.29: A control system with two inputs 


(a) Assuming that the initial conditions are zero, find the Laplace trans¬ 
form of the system output due to inputs: u\(t) = 2 e _t , u- 2 (t) = 4/?.(/), 
where h(t J is a unit step function. 

(b) Find y(t) by using the partial fraction expansion and the inverse 
Laplace transform. 

(c) Using MATLAB, find yi(t) = y(t) for u- 2 (t) = 0, and y 2 (t) = y(t) 
for ii i (/) = 0. Then, using superposition, find y(t) as a response to 
both u\(t) and u 2 (t). Plot the outputs, yi(t),y 2 (t), and y(t). 

2.6 Using the Laplace transform, find the transfer function of the electric 
network given in Figure 1.7, i.e. find V C3 {s)/Ei{s). 

2.7 Find the transfer functions Q(s)/V a (s) and Q,(s)/T^(s) of the armature- 
controlled DC motor given in Figure 2.6. 

2.8 Use formula (2.21) to find the closed-loop transfer function of a multivari¬ 
able system represented by 


r s+i 

s 

1 

■ r 

0 

r 

s 2 +6s+8 
— s 

s 2 +6s+8 

5 + 3 

II 

K 

S + l 

. 0 

Ls 2 +65+8 

5 2 +65+8 -1 

1 

s+2 . 


2.9 The block diagram of a control system is shown in Figure 2.30. 

(a) Reduce this system using block diagram algebra rules and find its 
transfer function Y(s)/U(s). 

(b) Draw the signal flow diagram of the system and find the transfer 
function by Mason’s rule. 

(c) Find the transfer function of the system using the MATLAB functions 
series and feedback repeatedly. 
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(d) Find the unit step response of this system using MATLAB. 



Figure 2.30: Block diagram of a control system 

2.10 Using block diagram algebra rules, simplify the block diagrams shown in 
Figure 2.31 and find their transfer functions (the matrix transfer function 
for multi-input multi-output case). 



Figure 2.31: Block diagrams 

(a) Repeat the procedure using Mason’s gain formula. 

(b) Find the transfer functions using MATLAB with 


G\(s ) 


3 

Kl+T)’ 


G 2 (s) 


5 

•s 2 + 2 


2.11 Verify that the block diagrams shown in Figure 2.32 arc equivalent, i.e. 
show that they have identical transfer functions Y(s)/U(s) and Y(s)/D(s). 
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Figure 2.32: Equivalent block diagrams 

2.12 Show that the transfer functions of the block diagrams given in Figures 
2.33a and 2.33b are identical, i.e. conclude that the block diagrams are 
equivalent. Note that this equivalence has been used to move from Figure 
2.14b to Figure 2.14c in Example 2.8. 



Figure 2.33: Equivalent block diagrams 

2.13 Show that both transfer functions Y(s)/Ui(s) and Y(s)/U- 2 (s) in the block 
diagrams given in Figures 2.34a and 2.34b are identical, i.e. conclude that 
the block diagrams are equivalent. Note that this equivalence has been used 
to move from Figure 2.14c to Figure 2.14d in Example 2.8. 




92 


TRANSFER FUNCTION APPROACH 



Figure 2.34: Equivalent block diagrams 

2.14 Find the transfer function between the nodes E 2 and E 5 , i.e. find the ratio 
E$(s)/E- 2 (s) for the signal flow graph given in Figure 2.17. 

2.15 SIMULINK can be used for drawing block diagrams. Once a block diagram 
has been drawn and a name to it given, we can obtain the state space 
model by invoking the SIMULINK function linmod ( ' block diagram 
name' ) . From the state space representation, the MATLAB function 
ss2tf produces the corresponding transfer function. Find the transfer 
functions for the block diagrams considered in Problems 2.9-2.13 by using 
the SIMULINK package. 

2.16 Omit the branch containing block G 5 in the block diagram presented in 
Figure 2.14a and find the system transfer function by using: 

(a) Mason’s rule. 

(b) Block diagram algebra. 

2.17 Using Mason’s gain formula, find the transfer functions of the systems 
whose signal flow graphs are shown in Figure 2.35. 

2.18 Consider the discrete-time linear system represented by the transfer function 

c 4 - 3c 3 + 5c 2 + 2c 

■ ■ ~~ 2U+5c 5 -3c 4 + c 3 —2c 2 + 3c— 1 

Use MATLAB in order to find: 

(a) The zeros and poles of the system. 
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(b) The inverse 2-transform of the transfer function. Hint: Use the 
residue function. 

(c) The closed-loop zeros and poles. 
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2.19 Find the inverse 2-transform for the following discrete transfer functions 

tt, ]= „ _ '2z 2 + s _ 

■■ z' 2 + 2z + 5’ ~ (~ + 0.4)(~ - 1)(~ - 0.2) 

Using MATLAB, find the first 10 samples of system outputs y(k J if the 
system input is u(k) = 1 ,k = 0,1,2,..., and u(k) = 0 for k < 0. 

2.20 For the digital control system shown in Figure 2.26 the plant transfer 
function is 

G(s) = —-— 

s(s + 1) 

(a) Find the discrete transfer function of the closed-loop system if the 
sampling interval is T = 0.5s, and the digital controller is 

Gd(z) = K p + y~^T 

(b) Using MATLAB, find the poles and zeros of the closed-loop system 
if K p = 1, Ki = 1. Plot the output response to a unit step input. 



